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ABSTRACT 

Array  detectors  based  on  correlator- 
- —type  structures  are  considered  for  robust 
sCaetection  of  weak  signals  with  noise  rep- 
^^resented  by  a generalized  Tukey-Huber 
QQcontami  nation  model.  The  performance  of 
__a  different  type  of  quadratic  limiter-sum 
C'JBetector  for  robust  detection,  obtained 
through  a consideration  of  the  locally 
^ ^optimum  detector  for  robust  detection, 
is  also  discussed. 

I . INTRODUCTION 

Robust  statistical  procedures  for 
signal  detection  problems  exhibit  per- 
formances which  are  fairly  stable  over 
classes  of  statistical  characterizations 
of  the  signals  and  noise.  They  are 
therefore  useful  in  implementing  detec- 
tors which  have  to  operate  under  some 
uncertainty  as  to  the  exact  nature  of 
the  signal  and  noise  processes.  The  . 
known  signal  in  additive  noise  problem 
has  been  considered  by  several  investi- 
gators (e.g.  [1]— [3]) . Techniques  for 
robust  detection  of  random  signals  from 
multiple  inputs  have  been  less  thorough- 
ly investigated,  although  nonparametric 
detectors  have  been  given  more  attention 
[4],  [5].  Robust  correlator-type  detec- 
tors have  been  considered  in  [ 6 ] and  [ 7 ] ; 
in  [6],  two- input  correlators  were  stud- 
ied and  some  limiter  structures  were  in- 
vestigated for  performance,  and  in  [7]  a 
general  robust  detection  problem  was  for- 
mulated. 

In  this  paper  we  show  that  if  atten- 
tion is  restricted  to  correlator-type 
structures  and  detector  performance  is 
characterized  in  terms  of  an  asymptotic 
measure,  the  efficacy  [8],  robust  detec- 
tors for  random  signals  in  additive  noise 
can  be  explicitly  defined.  We  will  use 
the  Tukey-Huber  contamination  noise  model 
to  define  the  class  of  noise  statistics; 
extensions  of  this  model  allow  us  to 
look  at  the  case  of  noise  processes  which 
are  correlated  between  the  several  input 
channels.  In  addition,  we  also  briefly 
^ — consider  for  robust  detection  a modified 
' ^“quadratic- sum  detector  which  performs 
CL.better  for  small  arrays. 

o 


II.  ROBUST  CORRELATORS 
A.  Preliminaries 

For  detecting  a random  signal  com- 
mon to  an  array  of  m channels,  we  assume 
that  a total  of  n vector  samples  X^, 

Xj,...,:^  are  taken,  where  X^  = (x^, 

Xi2,...,xim)  is  the  vector  of  outputs  at 

each  of  the  m channels  at  the  i-th  time 
instant;  thus  X^  is  the  output  of  the 

j-th  channel  at  the  i-th  sampling  time. 
If  0S^  represents  a random  signal  common 

to  the  array  at  the  i-th  sampling  time, 
we  can  describe  the  m-variate  probabil- 
ity density  function  fx | 0 of  X^  by 


fX  I 0 


/ fx(x  - 0s)dFg(s). 


In(l),Fg  is  the  distribution  function  of 

S^  and  is  assumed,  without  loss  of  general- 
ity, to  have  unit  variance;  the  vector  s 
is  an  m-component  vector  with  each  com- 
ponent equal  to  s;  and  fx  is  the  m-vari- 

ate  density  function  of  noise  across  the 
m channels.  With  this  notation,  f x | Q=f x- 

We  have  assumed  that  the  statistics  of 
the  signal  and  noise  are  stationary  in 
time,  and  we  will  also  assume  that  under 
the  null  hypothesis  HQ;  0=0,  the  time 

samples  X^  are  independent.  This  inde- 
pendence of  the  noise  samples  is  neces- 
sary for  the  rest  of  the  development. 

For  the  two-input  case  (m=2)  the 
general  correlator  detector  [ 6 ] is 
based  on  the  test  statistic 


S2  " iEiW<Xil>W(Xi2) ' ' (2) 

where  w is  a suitable  function.  Extend- 
ing this  for  the  m-input  case,  we  now 
consider  a correlator-array  test  statis- 
tic S_  defined  by 
m 

n m-1  m 

Sm  - I I I w(X. .)w(x.v)  . (3) 

m i-1  j-1  k=j+l  ^ lk 

We  will  consider  later  the  restriction  as 


, u;.e  . well  as  justification  for  the  assumption 
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sum  detector  for  robust  detection,  obtained  through  a consideration  of  the 
locally  optimum  detector  for  robust  detection,  is  also  discussed. 


of  ^his  structure  for  a multi-input  test 
statistic. 


2 - 


Under  rather  general  conditions, 
-the  ARE  of  two  detectors  can  be  ex- 
pressed as  a ratio  of  their  effica- 
cies [ 8 ] . Thus  our  measure  of  perform- 
ance will  be  the  efficacy  of  a correla- 
tor array.  The  efficacy  Effl  of  the  test 

statistic  S_  is  given  by 

JU 


*■- 


lim 

n-** 


n var.t  S } 
o m 


(4) 


where  e is  an  assumed  degree  of  contami- 
nation parameter  between  0 and  1.  The 
robust  detection  problem  is  then  to  find 
a "worst-case"  h,  say  hR,  and  corre- 
sponding density  f=f  R , for  which  the  op- 
timum detector  DR  is  derived.  DR  is 

then  the  most  robust  detector  if  for  any 
f,  the  performance  of  DR  (here  measured 

by  its  efficacy)  is  lower  bounded  by  its 
(optimum)  performance  for  fR 

He  now  consider  the  model  (8)  and 
sane  extensions  for  our  present  appli- 
cation. 


If  we  assume  now  that  fx,  the  m-variate 

density  of  the  array  noise  vector,  is  a 
product  of  m independent  and  identical 
symmetric,  univariate  densities  f,  i.e. 


f_(x)  - ir  f(x. ) , 
--  ~ i-1  1 


(5) 


we  have,  assuming  that  f (x)  is  suffi- 
ciently smooth  and,  (without  loss  of 
generality),  that  h(x_)  has  zero  mean 

under  HQ,  the  result  that 


i(m-l)  [/  w(x)f  (x)dx] 
""2  [/  w2(x)f (x)dx]2 


(6) 


On  the  other  hand,  if  w is  sufficiently 
Jth  and  bounded,  we  obtain 

,4 


t(m-l)  [/  w (x) f (x)dx] 
2 [/  w2(x)f(x)dx] 


7 • 


(7) 


Equation  (7)  is  the  result  which  we  will 
use  later. 

From  (6),  Schwarz's  inequality  im- 
mediately leads  to  the  conclusion  that 

w(x)  - c f ' (x)/f(x),  where  c is  any  con- 
stant, gives  the  optimum  correlator- 
array  when  f is  known.  We  now  consider 
robust  correlators  for  situations  where 
f,  and  more  generally  fx , is  not  com- 
pletely specified. 

The  contaminated  noise  model  ini- 
tially used  by  Tukey  and  Huber  for 
robustness  studies  [9],  [10]  (and  later 
used  to  formulate  signal  detection  prob- 
lems) provides  one  way  of  characterizing 
a class  of  densities  for  the  noise.  The 
univariate  models  which  have  been  used 
so  far  model  a density  f as  a mixture  of 
a known  nominal  component  g and  an  arbi- 
trary contamination  h; 

f(x)  -(l-e)g(x)  + e h(x)  , (8) 


B.  Robust  Correlators  for  Contaminated 
Noise 

The  first  noise  model  for  the  multi- 
input detection  problem  we  are  consider- 
ing is  obtained  by  combining  the  assump- 
tion of  (5)  with  the  model  (8) ; we  will 
call  the  result  Model  1.  The  resulting 
noise  class  contains  m-variate  densities 
for  the  array  noise  vector  of  the  form 


fx(x) 


it  [(l-e)g(x.)  + eh  (x- ) ] . (9) 
i-1  1 1 


We  will  assume,  for  Model  1 and  for  the 
other  models  we  will  consider,  that  the 
nominal  density  g is  symmetric,  -log  g(x) 
is  convex  (as  for  the  Gaussian  density) 
and  g' (x)  is  absolutely  continuous  (see 
[2  ]).  We  also  assume  for  Model  1 that 
h is  bounded  and  symmetric. 

Since  (9)  is  still  a product  of 

univariate  densities,  the  efficacy  E . 

m,i 

of  Sm  for  fx  a density  given  by  Model  1 


is 


m,  1 


<1 

m(m-l) 


-e)/w' (x)g(x)dx  + 

e/w  (x)h(x)dxj 


(10) 


-e)/wz(x)g(x)dx  + 


(1-E) 


e/w*(x)h(x)dx 


Comparison  of  this  expression  with  re- 
sults obtained  in  f 2 ] lead  us  directly 
to  the  following  theorem: 

Theorem  1 

For  the  class  of  noise  densities 
given  by  Model  1,  fR  defined  by 

l*|<a 


fR(x) 


m-e)g(x) 

lll-E)g(a) 


e)g (a) exp[-b ( |x|-a) 1 |x|>a 

(11) 

where  a,  b satisfy  /®g(x)dx  + 2g(a)/b  - 
(1-e)-1  and  -g'(a) /g (a)  « b,  and  the 


corresponding  nonlinearity 

wR(x)  - -fR(x)/fR(x)  (12) 

constitute  the  least  favorable  density 
and  the  maxi-min  correlator-array  non- 
linearity, respectively,  wR  giving  the 

most  robust  performance  in  terms  of 
efficacy. 

The  proof  is  quite  simple;  from  (11), 
we  see  that  hR  is  zero  in  (-a, a),  so  that 

arith  wR  replacing  w in  (10) , the  denomi- 
nator is  maximized  by  hR.  Also,  for  wR 
in  the  numerator,  hR  minimizes  the  value 
of  the  numerator,  wR (x)  being  zero  for 

|x|>a  and  non-negative  for  |x|<_a.  By 
(12) , we  conclude  that  wR  results  in  the 

optimum  correlator  for  fR,  and  therefore 

the  theorem  is  proved.  It  should  be 
noted  that  in  [ 7 ] the  nature  of  such  a 
saddle  point  solution  has  been  discussed 
for  a general  convex  family  of  noise  den- 
sities (here  generated  by  Model  1) . 

It  is  interesting  that  the  robust 
nonlinearity  wR  is  identical  to  the  ro- 
bust nonlinearity  for  the  known  signal 
in  additive  noise  detector  problem  con- 
sidered in  [2];  it  is  true,  however, 
that  for  the  known  signal  problem  the 
false-alarm  probability  was  kept  bounded, 
whereas  robustness  of  performance  in 
terms  of  efficacy  does  not  necessarily 
imply  a controlled  false-alarm  probabil- 
ity for  a fixed-threshold  test. 

We  now  consider  a different  exten- 
sion of  the  univariate  contamination 
model  of  (8) , which  does  not  give  chan- 
nel- to-channel  independence  for  the 
noise.  With  g and  h density  functions 
satisfying  the  same  conditions  as  in 
Model  1,  we  define  Model  2 to  give  in- 
var iate  noise  density  functions 

m m 

fY(x)  ■ (1-e)  n g(x. ) + e it  h(x. ).  (13) 
- i-1  1 i-1  1 

This  results  in  marginal  bivariate  densi- 
ties for  uncorrelated  random  variables; 
however,  it  is  clear  that  these  are  not 
independent.  It  is  easily  seen  that  the 
maximal  correlation  coefficient  [11], 
which  is  a much  more  suitable  measure  of 
dependence  between  random  variables,  may 
be  strictly  bounded  away  from  zero.  Thus 
for  Model  2,  e is  the  degree  of  contamina- 
tion and  dependence. 

For  Model  2,  it  follows  easily  that 


-the  efficacy  Em  2 of  the  correlator- 
array  statistic  (3)  is  given  by 

f(l-e) (/w1 (x) g (x)dx) 2 + 


_ m(m-l) 
Em-2  2 


c(/w1(x)h(x)dx)2 


(1-e) [/w2(x)g(x)dx]2  + 


e[/w2(x)h(x)dx]2  (14) 


A comparison  of  (14)  with  (10) , and  an 
examination  of  the  proof  of  Theorem  1 
allows  us  to  state  the  following  theorem: 

Theorem  2 

For  the  class  of  noise  densities 
given  by  Model  2,  wR  defined  in  (12)  re- 
sults in  a correlator-array  statistic 
which ‘is  robust  in  the  sense  that  its 

efficacy  E _ is  always  larger  than  or 
m,  z 

equal  to  the  efficacy  obtained  with  that 
h=hR  which  results  in  the  fR  of  (11) . 

The  proof  of  this  theorem  is  quite 
simple.  Note  that  we  cannot  conclude 
that  the  density  fx  R(x)  given  by 

m m 

fv  n(x)  - (1-e)  * gtx.)  + e n h (x.) 

and  constitute  a saddle-point  solution 

for  the  robust  detection  problem,  as  in 
Theorem  1.  This  is  because  it  has  not 
been  shown  that  wR  is  the  optimum  corre- 
lator-array nonlinearity  for  the  multi- 
variate noise  density  of  (14)  . The  re- 
sult is  nevertheless  very  useful  because 
it  gives  more  significance  to  the  robust 
nonlinearity  wR  of  Theorem  1;  with  wR, 

the  resulting  performance  remains  lower 
bounded  by  a computable  efficacy  for  all 
densities  in  Model  2 as  well.  Some  num- 
erical values  of  the  performances  for 
noise  densities  generated  by  Models  1 
and  2 are  shown  in  Table  I-III. 

The  numerical  values  shown  were  com- 
puted for  the  particular  case  where  g is 
a zero  mean,  unit  variance  Gaussian  den- 
sity. Table  I shows  the  lower  bounds  on 
the  efficacies  (normalized  by  dividing 

by  when  wR  is  the  correlator 

array  nonlinearity,  as  a function  of  the 
contamination  parameter  e.  The  two  lower 
bounds  (for  densities  in  the  two  models) 
are  very  close,  and  they  decrease  with 
increasing  e,  as  is  to  be  expected.  Of 
course,  the  lower  bounds  are  achieved 
for  a rather  special  kind  of  contasdna- 
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tion,  so  that,  for  example,  if  the  con- 
tamination density  is  also  unimodal, 
the  resulting  performance  will  usually 
be  considerably  better.  It  is  seen 
that  the  limiter-type  nonlinearity  can 
be  expected  to  give  reasonably  good  per- 
formance for  all  densities  generated  by 
Models  1 and  2,  for  small  contamination 
levels  e. 

In  Table  II,  normalized  efficacies 
are  shown  for  a specific  type  of  con- 
tamination in  Models  1 and  2;  the  den- 
sity function  h was  taken  to  be  Gaus- 
sian with  variance  a2,  with  a2>l,  to 
•observe  the  effect  of  broader  tailed 
contaminations . In  Table  III,  for  com- 
parison, the  same  contaminations  are 
applied  to  the  nominal  density  for  the 
correlator-array  which  is  optimal  for 
the  nominal  density,  i.e.  with  w(x)=x. 

It  is  clear  that  the  performance  of  the 
optimum  array  deteriorates  rapidly, 
especially  for  Model  2,  as  the  noise 
density  starts  getting  non-Gaussian. 

We  can  extend  Model  2 for  a con- 
taminated multivariate  nominal  density 
by  relaxing  the  requirement  that  the 
multivariate  contaminating  density 
m 

h(x)  is  a product  rr  h(x.)  of  univariate 
i-1  1 

densities.  Thus  we  consider  Model  3 for 
contaminated  multivariate  densities 
f x (x) , which  are  defined  for  this  model 

by  m 

f_(x)  - (1-e)  ir  f (x • ) + eh^x)  (15) 
- i-1  1 - ~ 

We  now  want  to  see  if  the  robust  corre- 
lator for  Models  1 and  2 also  has  a 
lower-bounded  performance  for  noise 
densities  in  this  larger  class  of  densi- 
ties than  that  generated  by  Model  2. 
Clearly,  some  restrictions  will  have  to 
be  imposed  on  h^. 

Suppose  we  assume  that  fx  and  hence 

hg  describes  spatially  homogeneous  and 

symmetrically  distributed  noise,  by 
which  we  mean  that  any  k-variate  margi- 
nal density  of  the  channel  noise  vector 
is  the  same  k-variate  function,  and  the 
value  of  the  noise  density  is  invariant 
under  changes  of  sign  for  any  k of  its 
m arguments,  for  k»l,2,...,m.  Under  the 
hosiogeneity  and  symmetry  assumptions,  we 
can  show  that  the  robustness  of  the  cor- 
relator-array nonlinearity  wR  for  Model 

2 also  holds  for  the  larger  class  of 
noise  densities  of  Model  3.  In  spite  of 
the  conditions  of  symmetry  and  homo- 
geneity, Model  3 generates  a strictly 


larger  class  of  densities  than  does 
Model  2;  for  example  one  can  include 

in  Model  3 which  are  themselves  mixtures 

m 

of  contaminations  of  the  form  it  h(x.). 

i-1  1 

Finally,  it  should  be  noted  that  it 
is  also  possible  to  show  that  a lower 
bound  on  performance  exists  for  the  cor- 
relator-array using  w even  without  the 

K 

spatial  homogeneity  assumption.  This  is 
a direct  extension  of  the  above  results. 
For  example,  in  a geometrically  linear 
array,  it  is  more  reasonable  to  assume 
that  the  bivariate  inter-channel  noise 
densities  could  depend  on  the  spacing 
between  the  channel  receivers.  Such 
consideration*  will  modify  the  expres- 
sion for  the  derivative  of  the  mean  and 

variance  of  the  test  statistic,  the 

terms  now  not  being  all  identical.  It 
will  still  be  true,  however,  that  the 
performance,  in  terms  of  efficacy,  of 
the  correlator-array  using  wR  is  lower 

bounded  as  before. 


III.  ROBUST  ARRAY  FROM 
LOCALLY-OPTIMUM  ARRAYS 

All  the  previous  results  were  ob- 
tained for  the  case  where  the  signal 
detection  structure  was  constrained  to 
be  of  the  correlator  type.  One  of  the 
consequences  of  this  is  that  no  results 
were  obtained  for  the  one-channel  random 
signal  problem.  Theorem  1 tells  us  that 
wR  is  the  most  robust  nonlinearity  sub- 
ject to  this  correlator  structure  con- 
straint. We  will  now  consider  briefly 
a robust  detector  structure  motiviated 
by  the  form  of  the  locally -optimum  array 
detector. 


Applying  the  generalized  Neyman- 
Pearson  lemma  [8],  we  find  that  the 
locally-optimum  detector  for  a random 
signal  in  additive  noise  in  m channels 
is  based  on  the  test  statistic  T given 
by 


nr  m f " (X.  .) 

z h r kx-t 

i*l  l j*l  tlxij' 


“I1  ® f (xij}  f (Xik)  1 

j*l  k-j+1  f(Xij*  f {^ik*  J 
«• 


(16) 


assuming  f exists  almost  everywhere. 
It  is  clear  from  (16)  that  the  optimum 
array  detector  for  weak  signals  is  not 
a correlator  in  general;  in  fact,  T 


reduces  to  a correlator-array  only  for 
the  double-exponential  density,  giving 
the  polarity-coincidence  array  [4],  If 
m is  fairly  large,  the  m second  deriva- 
tive terms  in  (16)  may  be  ignored,  be- 
ing a minor  contribution  in  comparison 
with  the  m(m-l)  correlation  terms. 

'Thus  for  large  m,  a correlator-struc- 
ture assumption  is  quite  reasonable. 

The  correlator-array  detector  has  gen- 
erally received  considerable  attention 
because  of  its  relative  ease  of  imple- 
mentation, although  it  is  known  that  in 
the  two-input  case  [12],  for  example, 
its  ARE  relative  to  the  locally-optimum 
detector  is  quite  poor.  This  is  corre- 
spondingly true  for  robust  detection 
based  on  the  correlator-structure  as 
opposed  to  the  locally-optimum  struc- 
ture. 


Consider  the  case  of  Gaussian  noise 
with  unit  variance;  then  the  statistic  T 
becomes,  after  dropping  constant  terms. 


n m 

E ( E X.  .) 

i«l  tel  13 

•<c*  ■ 


n 

- I 

i-1 


m , 
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This  can  also  be  written  as 
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where  f here  is  the  unit-variance  Gaus- 
sian density.  The  structure  of  (18) 
leads  us  to  investigate  the  robustness 
properties  of  the  test  statistic  TR 
defined  by 
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not  very  different.  The  TR  of  (19)  was 

motivated  purely  by  the  previous  results, 
and  it  is  not  the  most  robust  statistic 
in  that  a saddlepoint  for  performance  is 
not  obtained  with  TR. 

IV.  SUMMARY 

Robust  correlator  structures  were 
considered  for  random  signal  detection 
in  additive  noise.  The  criterion  of 
performance  used  was  the  efficacy  of  a 
detector  test  statistic,  a measure  of 
the  weak-signal  detection  performance 
of  the  statistic.  For  a direct  multi- 
variate extension  of  the  Tukey-Huber 
contaminated  noise  model,  it  was  shown 
that  the  max-min  robust  nonlinearity  is 
the  same  as  the  robust  nonlinearity  for 
known  signals  in  contaminated  additive 
noise.  The  contamination  noise  model 
was  then  generalized  to  allow  non-inde- 
pendent multivariate  contaminations.  It 
was  shown  that  the  same  robust  nonlin- 
earity provides  stable  performance  for 
the  resulting  class  of  noise  densities. 
Finally,  we  considered  the  performance 
of  an  array  detector  for  robust  detec- 
tion for  small  arrays.  The  quadratic 
limiter-sum  detector  for  Gaussian  noise 
was  found  to  have  better  stable  perform- 
ance for  contaminated  noise  in  small 
arrays,  and  its  performance  was  close 
to  that  of  the  robust  correlators  in 
moderate-sized  arrays. 
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1 

4 

9 
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0.787 
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Efficacy, 
Model  2 

0.913 
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